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Abstract 

The associated production of a Higgs boson with a b quark is a discovery channel for the hghtest 
MSSM neutral Higgs boson. We consider the SUSY QCD contributions from squarks and gluinos 
and discuss the decoupling properties of these effects. A detailed comparison of our exact 0{as) 
results with those of a widely used effective Lagrangian approach, the A;, approximation, is pre- 
sented. The Afe approximation is shown to accurately reproduce the exact one-loop SQCD result 
to within a few percent over a wide range of parameter space. 
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I. INTRODUCTION 



Once a light Higgs-like particle is discovered it will be critical to determine if it is the 
Higgs Boson predicted by the Standard Model. The minimal supersymmetric Standard 
Model (MSSM) presents a comparison framework in which to examine the properties of 
a putative Higgs candidate. The MSSM Higgs sector contains 5 Higgs bosons-2 neutral 
bosons, h and if, a pseudoscalar boson. A, and 2 charged bosons, H"^. At the tree level 
the theory is described by just 2 parameters, which are conveniently chosen to be Ma-, the 
mass of the pseudoscalar boson, and tan/3, the ratio of vacuum expectation values of the 
2 neutral Higgs bosons. Even when radiative corrections are included, the theory is highly 
predictive [iHS I . 

In the MSSM, the production mechanisms for the Higgs bosons can be significantly 
different from in the Standard Model. For large values of tan/3, the heavier Higgs bosons, 
A and if, are predominantly produced in association with h quarks. Even for tan /3 ~ 5, 
the production rate in association with h quarks is similar to that from gluon fusion for 
A and H pro duct ion For the lighter Higgs boson, h, for tan/3 > 7 the dominant 
production mechanism at both the Tevatron and the LHC is production with h quarks for 
light Ma (^200 GeV), where the hbh coupling is enhanced . Both the Tevatronjj] and the 
LHC experimentsy] have presented limits Higgs production in association with h quarks, 
searching for the decays h — t^t~ and Ur . These limits are obtained in the context of the 
MSSM are sensitive to the 6-squark and gluino loop corrections which we consider here. 

The rates for hh associated production at the LHC and the Tevatron have been extensively 



studied [8 
flavor number parton schemes 



- |l8[ | and the NLO QCD correction are well understood, both in the 4- and 5- 

9, 111, [l^. In the 4- flavor number scheme, the lowest order 
processes for producing a Higgs boson and a h quark are gg hbh and qq — 66/ij8l, Q, IitI . 
In the 5— flavor number scheme, the lowest order process is hg — hh (bg — )■ bh). The 
two schemes represent different orderings of perturbation theory and calculations in the two 
schemes produce rates which are in qualitative agreement UJ, ll9| . In this paper, we use 



the 5-flavor number scheme 
electroweak corrections 



21 



'or simplicity. The resummation of threshold logarithms 20|] , 



22| and SUSY QCD corrections 



23| have also been computed for 



^ The expected sensitivities of ATLAS and CMS to b Higgs associated production are described in Refs. 
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bh production in the 5— flavor number scheme. 

Here, we focus on the role of squark and gluino loops. The 
QCD corrections to the bbh vertex, both for the decay h ^ bb 



24-: 



properties of the SUSY 



bb->h 
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27| and the production, 



27H29|. were computed long ago. The contributions from b squarks and gluinos to 



the lightest MSSM Higgs boson mass are known at 2-loops[30|, ISll, while the 2-loop SQCD 
contributions to the bbh vertex is known in the limit in which the Higgs mass is much smaller 



than the squark and gluino masses 



32 



33|. The contributions of squarks and gluinos to the 



on-shell bbh vertex are non-decoupling for heavy squark and gluino masses and decoupling 
is only achieved when the pseudoscalar mass. Ma, also becomes large. 

An effective Lagrangian approach, the A;, approximation 25|, l26| , can be used to approxi- 
mate the SQCD contributions to the on-shell bbh vertex and to resum the (a^ tan /3/Msusy)"' 
enhanced terms. The numerical accuracy of the Af, effective Lagrangian approach has been 
examined for a number of cases. The 2— loop contributions to the lightest MSSM Higgs 
boson mass of 0{ab<ys) were computed in Refs. [30I and 31|, and it was found that the 
majority of these corrections could be absorbed into a 1— loop contribution by deflning an 
effective b quark mass using the A;, approach. The sub-leading contributions to the Higgs 
boson mass (those not absorbed into A5) are then of 0{1 GeV). The A^ approach also 
yields an excellent approximation to the SQCD corrections for the decay process h — > 66 j^^]. 
It is particularly interesting to study the accuracy of the A;, approximation for production 
processes where one of the b quarks is off-shell. The SQCD contributions from squarks and 
gluinos to the inclusive Higgs production rate in association with b quarks has been studied 
extensively in the 4FNS in Ref. js^], where the the lowest order contribution is gg — ?■ bbh. 
In the 4FNS, the inclusive cross section including the exact 1-loop SQCD corrections is 
reproduced to within a few percent using the A;, approximation. However, the accuracy of 
the Ab approximation for the MSSM neutral Higgs boson production in the 5FNS has been 
studied for only a small set of MSSM parameters in Ref. |23|. The major new result of this 
paper is a detailed study of the accuracy of the Af, approach in the 5FNS for the bg — > bh 
production process. In this case, one of the b quarks is off-shell and there are contributions 
which are not contained in the effective Lagrangian approach. 

The plan of the paper is as follows: Section 2 contains a brief review of the MSSM Higgs 
and b squark sectors and also a review of the effective Lagrangian approximation. The 
calculation of Ref. j23l] is summarized in Section 2. We include SQCD contributions to 
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bh production which are enhanced by rrib tan /3 which were omitted in Ref . 23| . Analytic 
results for the SQCD corrections to bg — ?■ bh in the extreme mixing scenarios in the b squark 
sector are presented in Section 3. Section 4 contains numerical results for the a/s = 7 TeV 
LHC. Finally, our conclusions are summarized in Section 5. Detailed analytic results are 
relegated to a series of appendices. 



II. BASICS 



A. MSSM Framework 



In the simplest version of the MSSM there are two Higgs doublets, and H^, which 
break the electroweak symmetry and give masses to the W and Z gauge bosons. The neutral 
Higgs boson masses are given at tree level by, 

1 



MI + MIt \Km\ + M|)2 - AMlMl cos2 2/3 



and the angle, a, which diagonalizes the neutral Higgs mass is 

'Ml + Ml' 



tan 2a = tan 2/3 



(1) 



(2) 



In practice, the relations of Eqs. [T] and [2] receive large radiative corrections which must be 
taken into account in numerical studies. We use the program FeynHiggs 34 to generate 
the Higgs masses and an effective mixing angle, Oe//, which incorporates higher order effects. 

The scalar partners of the left- and right- handed b quarks, bi and b^, are not mass 
eigenstates, but mix according to. 




(3) 



The b squark mass matrix is. 




(4) 



and we define, 



Xfc = Ah — fi tan /3 



m 



R 



MQ + ml + Ml cos 2l3{ll - Qb sin^ Ow) 
Ml + ml + Ml cos 2/3gbsin2 Ow 



(5) 



Mq^d are the soft SUSY breaking masses, /| 



-1/2, and Qh 



-1/3. The parameter 



is the trihnear scalar couphng of the soft supersymmetry breaking Lagrangian and /i is 
the Higgsino mass parameter. The b squark mass eigenstates are bi and 62 and define the 
6-squark mixing angle, 6^ 

bi = cos 9bbL + sin Ofybfi 
62 = — sin Obbi + cos 9bbR . 

(6) 

At tree level, 

. ~ 2mfe(v4b - /itan/3) 

Mf-M^ (7) 

and the sbottom mass eigenstates are, 

' ' (8) 



M? - 

bi,b2 



B. Ab Approximation: The Effective Lagrangian Approach 



Loop corrections which are enhanced by powers of as tan /3 can be included in an effective 
Lagrangian approach. At tree level, there is no tp^bj^Hu cou plin g in the MSSM, but such a 



coupling arises at one loop and gives an effective interaction 

A, 



25 



27| 



'eff 



-Xb^L Hd + 



* " ' tan /3 

Eq. [9] shifts the b quark mass from its tree level value, ^ 



Hu bR + h.c. 



nib —^{^ 



V2 



(9) 



(10) 



and also implies that the Yukawa couplings of the Higgs bosons to the b quark are shifted 
from the tree level predictions. This shift of the Yukawa couplings can be included with an 
effective Lagrangian approach (26|, |27|, 

rrib ( 1 \ ( sina\ ( ^ Ab 



The 

tan/3 



L,fj = -^^-^ 1 ^ — \bbh. (11) 

vsM\l + AbJ\ COS (3 J \ tan/3tanay ^ ^ 

lagrangian of Eq. [11] has been shown to sum all terms of (9(astan/3)" for large 



25 



26|.^ This effective Lagrangian has been used to compute the SQCD corrections 



^ The neutral components of the Higgs bosons receive vacuum expectation values: (ff") = (-ff-S) — 
^ vsM = (^/2G'F)"^/^ vi = wsM cos /3 

^ It is also possible to sum the contributions which are proportional to Af,, but these terms are less important 
numerically 27 1 . 
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to both the inclusive production process, bb h, and the decay process, h — >■ 65, and yields 
results which are within a few percent of the exact one-loop SQCD calculations 



27 



37|. 



The expression for Ah is found in the limit rrih << Mh„ Mz << Af^^, M^^, Mg . The 1-loop 
contribution to from sbottom/gluino loops is[25|, \2m, l38| 



A, 



where the function J (a, b, c) is. 



-M5/itan/3J(M,;,M,-^,M5) 



(12) 



/(a, 6, c) 



[a^ - b-) ib- - c^) (a^ - c^) 1 ) + ( ^ ' + 



(13) 

and as{ns) should be evaluated at a typical squark or gluino mass. The 2— loop QCD 



corrections to Ab have been computed and demonstrate that the appropriate sca^ 



e at which 



to evaluate Af, is indeed of the order of the heavy squark and gluino masses |32l. |33|. The 
renormalization scale dependence of A;, is minimal around /io/3, where /iq = {Mg + m^^ + 
'mi^)/3. In our language this is a high scale, of order the heavy SUSY particle masses. 
The squarks and gluinos are integrated out of the theory at this high scale and their effects 
contained in A;,. The effective Lagrangian is then used to calculate light Higgs production 
at a low scale, which is typically the electroweak scale, ~ 100 GeV. 

Using the effective Lagrangian of Eq. [9l which we term the Improved Born Approximation 
(or Afy approximation), the cross section is written in terms of the effective coupling, 

1 \ Afe 



9mi = 9bbh 



where 



dbbh — — 



1 + Ab J \ tan (3 tan a 
sin a\mb{fJ.R) 



cos /3 J VSM 



(14) 



(15) 



We evaluate WihifJ'R) using the 2— loop MS value at a scale fiR of 0{Mh), and use the value of 
determined from FeynHiggs. The Improved Born Approximation consists of rescaling 
the tree level cross section, o"o, by the coupling of Eq. [1] 



O'IBA 



9bbh 
Qbbh 



CTq. 



(16) 



The Improved Born Approxirnation has been shown to accurately reproduce the full SQCD 
calculation of pp — )■ tbH^ 



39, 4o|- 



This is the approximation used in Ref. [19i] to include the SQCD corrections. 
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FIG. 1: Feynman diagrams for g{qi) + h[q2) — )■ + h{ph)- 



The one-loop result including the SQCD corrections for hg — )■ hh can be written as, 



CSQCD = 0-IBa(^^ + ^SQCD^ , 



(17) 



where Asqcd is found from the exact SQCD calculation summarized in Appendix B. 

The Improved Born Approximation involves making the replacement in the tree level 
Lagrangian, 



rrib 7- . (18) 

Consistency requires that this substitution also be made in the squark mass matrix of Eq. 



^41 



42| 







V 



m 



1+At 



X, 



™2 



m 



R 



J 



(19) 



The effects of the substitution of Eq. [TH] in the 6-squark mass matrix are numerically 
important, although they generate contributions which are formally higher order in as- Eq s. 
and [19] can be solved iteratively for M^^, M^^ and A;, using the proceedure of Ref. 41p. 



C. SQCD Contributions to gb hh 



The contributions from squark and gluino loops to the gh — )■ hh process have been com- 



puted in Ref. 



23| in the mh = limit. We extend that calculation by including terms which 



are enhanced by ruf, tan /3 and provide analytic results in several useful limits. 

The tree level diagrams for g{qi) + b{q2) — > b{pb) + h{ph) are shown in Fig. [H We define 



We use FeynHiggs only for calculating Mh and sinae//. 
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the following dimensionless spinor products 

^ ujPb) + (92) 

* s 

* t 



= fb 



u 
rrib 

u{pb) iiu{q2) 



rribt 



Mr = g,-^iMil!iM, (20) 
nibS 

where s = {qi + q2Y,t = {pb — qiY and u = {pf, — ga)^- In the = limit, the tree level 
amplitude depends only on Mjf and Mj^, and Mf is generated at one-loop. When the effects 
of the b mass are included, M^*, Mg*, and M4 are also generated. 
The tree level amplitude is 

A^^p lo = -9s9bbk {T^ap e,{q,) {M^ + , (21) 



= -^^9s9bbH in^, ^^K'M ■ (22) 



and the one loop contribution can be written as 

Air 

In the calculations to follow, only the non-zero Xj coefficients are listed and we neglect 
terms of 0{ml/s) if they are not enhanced by tan/3. 

The renormalization of the squark andgluino contributions is performed in the on-shell 



scheme and has been described in Refs. 



23 



32 



43J. The bottom quark self-energy is 
ip) = ^(^^Up") - ^tip'H^ + mb^iip') ■ (23) 

The b quark fields are renormalized as 6 — > a/ b and = a/1 + SZ^ . The contribution 
from the counter-terms to the self-energy is, 

Sr(p) = S,(j9) + (5S,(p) 

(p) = i [bZX - bZti^) - m^bZX - (5m, . (24) 

Neglecting the 75 contribution, the renormalized self-energy is then given by 

sr (p) = {i- ^b) i^Up') + szr) + m, f sf (p^) + Er (p^) - ^) . (25) 



The on-shell renormalization condition implies 



lim 



The mass and wavefunction counter-terms are^ 
6mb 




0. 



(26) 
(27) 



nib 



1=1 



Ml 



-ly '-^s^-,Bo - B, 



0;MlMl 



(28) 



d 



El + 2771^5; - {-lY2mbMgS^iB'Q 



O 9 



(29) 



where we consistently neglect the h quark mass if it is not enhanced by tan /3. The Passarino- 
Veltman functions Bq ^0; M|, M-^ j and Bi ^0; M|, M| j are defined in Appendix A. Using 
the tree level relationship of Eq. [71 the mass counterterm can be written as, 
5mb 2as{^ii 

' ^ i=l 



3, -/,A,/(M,^, M,^, Mg) -Ab- ^ J2 (O' ^l) ■ (30) 
The external gluon is renormalized as — )■ \fZlg^ = a/1 + SZ^g^ and the strong cou- 



pling renormalization is gs — )■ Zggs with 5Zg = — 5Z3/2. We renormalize gs using the MS* 
scheme with the heavy squark and gluino contributions subtracted at zero momentum 4J], 



SZ. 



Qs(/ii 

4:71 



6 ''V ^< 



:r(l + e) 



(31) 



In order to avoid overcounting the effects which are contained in g^^^ to 0{as), we need 
the additional counterterm, 

1 



Ab 1 



tan (3 tan a 



The total contribution of the counterterms is. 



acT — C^IBA 



(2SZ, 



6Z. + 26Z, + 2 



6mh 



+ 26, 



2a 



szl 



6mh 



nib 



+ 5, 



CT 



(32) 



(33) 



The tan /3 enhanced contributions from A;, cancel between Eqs. [30] and [321 The expressions 
for the contributions to the Xi, as defined in Eq. [221 are given in Appendix B for arbitrary 
squark and gluino masses, and separately for each 1— loop diagram. 



'26 
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III. RESULTS FOR MAXIMAL AND MINIMAL MIXING IN THE &-SQUARK 
SECTOR 



A. Maximal Mixing 

The squark and gluino contributions to bg — )■ bh can be examined analytically in several 
scenarios. In the first scenario, 

We expand in powers of ^iXt ■ ^'^^^ ^^^^ sbottom masses are nearly degenerate, 



^ 2 



Ml +Mi 

bi 62 



- = " ' 1 + ^-^^ % 1. ^' «Ml. 35) 



This scenario is termed maximal mixing since 

sin 2Bb^\- ^ ^ / ^ • 36 

We expand the contributions of the exact one-loop SQCD calculation given in Appendix B 

1 and assuming Ms ~ Mg ~ /i ~ ~ niL ~ 
m/j >> Mw, Mz, Mh ~ Mew- In the expansions, we assume the large tan /3 limit and take 
mfttan /3 ~ 0{Mew)- This expansion has been studied in detail for the decay h — )■ bb, with 
particular emphasis on the decoupling properties of the results as Ms and Mg — t- ooj28|. 
The SQCD contributions to the decay, h — > bb, extracted from our results are in agreement 



with those of Refs. 



28 



The final result for maximal mixing, summing all contributions, is, 

= -9sT gbbhMN 1 + + — -JT2^l^r. 

I \9bbhJ^ax 47r 

'Sgbbh 



-gsT^'gbbhMN 1 



9bbh 



A^ ^ -9sT^gbbHM^[l + ^^Xl^ 



-9sT^9bbHM,^ (-^^) • (37) 
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The contribution which is a rescahng of the bbh vertex is, 



5g, 



bbh 



5g, 



bbh 



(1) 



+ 



9bbh / rnax \ 9bbh / ^nax 

where the leading order term in Mew l^s is 



Ibbh 



Qbbh 



(2) 



(38) 



bbh 



(1) 



37T 



Ml 



(39) 



9bbh / max 

with Yb = Ab + fi cot a and R = MJMs- Eq. |39]only decouples for large Ms if the additional 



limit Ma — > oo is also taken 



23 



28|. In this limit, 



2/iM| 

"mT 



tan /3 COS 2/3 + C 



M\ 



(40) 



The subleading terms of 0{M'^t^ / M^) are. 



^9. 



bbh 



9bbh 



(2) 



R 



Stt 



M-gYb 



Ml 



Ml 1 
-72J3 i^) 



12M| 



2(1 + A,)2M4 



nilXbYb 



2(1 + A;,)2M| 



/3"'(^) 

3/,(i?)+(^-l)/.(i?) 



(41) 



The functions fi{R) are defined in Appendix C. 
The ' 

obtained from the effective Lagrangian. We find 



^■i 1 ^ terms in Eq. |37] are not a rescaling of the lowest order vertex and cannot be 



8Kr. 



~R- 



Yb_ 

■2Mc 



+ \h{R) 



(42) 



The Suraax term is 0(V) in Mew /Ms and has its largest values for small R and large ratios 
of Yb/Ms, as can be seen in Fig. [2l Large effects can be obtained for Y^/Ms ~ 10 and 
Mg « Ms- However, the parameters must be carefully tuned so that A^/Ms ^ 1 in order 



not to break color 



45|. 



The amplitude squared, summing over final state spins and colors and averaging over 
initial state spins and colors, including one-loop SQCD corrections is 



o 9bbh 



u^ + M^ 



St 



59 



bbh 



Qbbh 



a^Ml 
27r M| " 



.(43) 



We use the shorthand, cp = cos/3, Sa+p = sin(a + etc. 
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FIG. 2: Contribution of 5 

f^max defined in Ec[. 1421 as a- function of i? = Mg/Ms- 



Note that in the cross section, the 5Kmax term is not enhanced by a power of s and gives a 



contribution of O ( 



Expanding in the maximal mixing hmit 
By comparison with Eq. [HI 



^^tan/3A(i?) + 0^^^^ 



(44) 



\A\ 



St 



^9^ 



bbh 



Qbbh 



(2) 



2tt 



Ml 



o 



M 



EW 



, a: 



(45) 



Note that the mis-match in the arguments of as in Eqs. HH and |45] is higher order in as 

(2) 

than the terms considered here. The {6 g /g bbh) max and SK^ax terms both correspond to 
contributions which are not present in the effective Lagrangian approach. These terms are, 
however, suppressed by powers of M^^r/Ml and the non-decouphng effects discussed in 



Refs. 



28 1 and 27 1 are completely contained in the g^'l^ term. 
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B. Minimal Mixing in the b Squark Sector 



The minimal mixing scenario is characterized by a mass sphtting between the b squarks 

' -M? |~ Ml 

1 b2 ' ^ 



which is of order the b squark mass, | — M? |~ M|. In this case 



(46) 



and the mixing angle in the b squark sector is close to zero, 



62 

The non-zero subamplitudes are 



As = -9sT^gbbHM^{ + ^^^^5/.. 

A = -gsT%bHMt{l + (^-^ 

I V 9bbh 

^1 = -gsT^QbbHM,' (-^^^] ^^^^n . (48) 

V 27rM^ / 

Expanding the exact one-loop results of Appendix B in the minimal mixing scenario. 



^f~\R,) + f,iR,) 



h ^^^^ [Sh{R,, R,, 1) + R,, 2) j , 

(49) 

where Ri = Mg/M^, and the functions fi{Ri) and hi{Ri, R2,n) are defined in Appendix C. 
The Sumin function is shown in Fig. El For large values of Y^/Mg it can be significantly 
larger than 1. 

As in the previous section, the spin and color averaged amplitude-squared is. 



-J ,2 _ £VA-K7^/ 2 ^ 

^ I mm — o yybbh) 



St 



1 + 2 



Jbbh 



9bbh 



27r ^- — Ml 



^gbbh\ ( ^9bbh\^^^ _|_ ( Sgbbh \ ^"^^ ^g-j^^ 



with. 



9bbh / min V 9bbh / min \ 9bbh / 

The leading order term in Mew /Ms is 

(1) 



(- — = — -B2 ^hi{Ri,R2,0) . [52) 

\ 9bbh ) <^7r Mg R^ - R2 
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FIG. 3: Contribution of SKmin defined in Eq. [39] as a function of Ri = Mg/Mj^ . 



The subleading terms are O 



M'i 



bbh 



Qbbh 



(2) 



An 1 3 AM^ 

^ bi2 



/l2 {Ri,R2)Ml 



+ 



AM^ 

012 



25 



AM^ 

012 



3 



5 



3M? 





AM? 

012 



O Srv 



A 



3/i (i?) - /2 (i?) 



3M^ 





'am? 

012 



5 



M? 

b 



AM? 

012 



8 mlX.Y, 
'2, AM? (1 + At 

012 



3/i (i?) - /2 {R) 



3M? 

b 



The symmetric and anti-symmetric functions are defined, 



Sif{R,M-,) ^ - 



AifiR,Mj, 



f{R^,M-,J + f{R,,M-,;, 
fiRuM-J-fiR2,Mi;^ 



3/i (R) - /2 (R) \ 2MgX,^ ( f\ (R) 



M? 
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and AM^ 

012 



— M~ . The remaining functions are defined in Appendix C. 



By expanding Aj, in tlie minimal mixing limit, we find the analogous result to that of the 
maximal mixing case, 



A 



Q yybbh) 



St 



1 + 2 



^9- 



bbh 



Qbbh 



(2) 



27r ™"M| ^ 



EW 



(55) 



The contributions which are not contained in ajBA are again found to be suppressed by 

Mew 



o 



Ms 



IV. NUMERICAL RESULTS 



We present results for pp h{h)h at = 7 TeV with pxb > 20 GeV and | 77^ |< 
2.0. We use FeynHiggs to generate Mh and sinoe// and then iteratively solve for the 



b squark masses and A;, from Eqs. [I2] and [191 We evaluate the 2-loop MS b mass at 
fiji = Mh/2, which we also take to be the renormalization and factorization scales^. Finally, 
Figs m O O and [7] use the CTEQ6m NLO parton distribution functionsj46j. Figs. HI [5] 
and [6] show the percentage deviation of the complete one-loop SQCD calculation from the 
Improved Born Approximation of Eq. [16] for tan /3 = 40 and tan /3 = 20 and representative 
values of the MSSM parameters'^. In both extremes of b squark mixing, the Improved 
Born Approximation approximation is within a few percent of the complete one-loop SQCD 
calculation and so is a reliable prediction for the rate. This is true for both large and small 
Ma- In addition, the large M5 expansion accurately reproduces the full SQCD one-loop 
result to within a few percent. These results are expected from the expansions of Eqs. [45] 
and [551 since the terms which differ between the Improved Born Approximation and the 
one- loop calculation are suppressed in the large M5 limit. 

Fig. [3 compares the total SQCD rate for maximal and minimal mixing, which bracket 
the allowed mixing possibilities. For large Ms, the effect of the mixing is quite small, while 
for Ms ~ 800 GeV, the mixing effects are at most a few fb. The accuracy of the Improved 
Born Approximation as a function of m/j is shown in Fig. [HI for fixed M^,^, and rriL. As 



^ A;, is evaluated using as{Ms)- 

10 pigg m [5] and [51 do not include the pure QCD NLO corrections 17] . 
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FIG. 4: Percentage difference between the Improved Born Approximation and the exact one-loop 
SQCD calculation oipp — )• hh for maximal mixing in the 6-squark sector at \fs = 7 TeV, tan /3 = 40, 
and Ma = 1 TeV. 



rriR is increased, the effects become very tiny. Even for light gluino masses, the Improved 
Born Approximation reproduces the exact SQCD result to within a fev^^ percent. 

In Fig. [9l we show the scale dependence for the total rate, including NLO QCD and SQCD 
corrections (dotted lines) for a representative set of MSSM parameters at ^/s = 7 TeV. The 
NLO scale dependence is quite small when fiR = fip ^ Mh. However, there is a roughly ~ 5% 
difference between the predictions found using the CTEQ6m PDFs and the MSTW2008 NLO 



PDFs 



47l |. In Fig. [TOl we show the scale dependence for small fip (as preferred by 



16|), 



and see that it is significantly larger than in Fig. [91 This is consistent with the results of 



19|, 



29|. 



V. CONCLUSION 



Our major results are the analytic expressions for the SQCD corrections to b Higgs 
associated production in the minimal (Eqs. |4T], |42] and l45l) and maximal (Eqs. |49l [53l 
and ES]) b squark mixing scenarios for large tan/3 and squark masses, M^. These results 
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FIG. 5: Percentage difference between the Improved Born Approximation and the exact one-loop 
SQCD calculation oipp — t- bh for maximal mixing in the 6-squark sector at ^/s = 7 TeV, tan /3 = 20, 
and Ma = 250 GeV. 



clearly demonstrate that deviations from the Ab approximation are suppressed by powers 
of {Mew /Ms) in the large tan/3 region. The A;, approximation hence yields an accurate 
prediction in the 5 flavor number scheme for the cross section for squark and gluino masses 
at the TeV scale. As a by-product of our calculation, we update the predictions for b Higgs 
production at = 7 TeV. 
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Appendix A: Passarino-Veltman Functions 

The scalar integrals are defined as: 



2^ 



167r2 



(27r)" N0N1N2NS 



where, 



N2 = ik + pi+p2f-M? 



2 

Ns = {k + p^+p2+p^Y-Ml. (57) 



d"A; 1 



dJ'k 1 



Do{pIpIpIpI {P1+P2Y, {p2+P2f\MlMlMlMl) 



(56) 



18 



7 TeV LHC, bg^bh 



50 



tan B=40, Li=M =lTeV, M =2M , = 






-\ 






\ 

\ 


Maximal Mixing 




- \ 

\ 


Minimal Mixing 




_ \ 
I \ 
\ V 
\ \ 
\ \ 
\ \ 
\ \ 
\ \ 

\ \ 

\ ^ 
\ \ 
\ \ 
\ \ 
\ N 

\ N 
\^ S 




- 









600 800 1000 1200 1400 1600 1800 2000 

Mj, (GeV) 



FIG. 7: Comparison between the exact one-loop SQCD calculation for pp — t- bh for minimal and 
maximal mixing in the b squark sector at = 7 TeV and tan /3 = 40. The minimal mixing curve 
has rriji = y/2Ms and 06 ~ 0, while the maximal mixing curve has m^j = Ms and 0;, ~ f . 



The tensor integrals encountered are expanded in terms of the external momenta Pi and 
the metric tensor g^'^ . For the two-point function we write: 



167r2 J {2tiY NqN: 



' ptBM.MlMl), (58) 



IGtt 

while for the three-point functions we have both rank-one and rank-two tensor integrals 
which we expand as: 

C^{pIpI{Vi+V2?;MImImI) = p'^C^.+p^.Cu, 
C^''{pl,pl{pi+P2?;MlM!,Mi) = pIp'.C^i + pi^p',C22 

+ W2+P\p'2)C2, + rC2,, (59) 

where: 

^C^C-n(pf.pI.te+P.)-A«.Mf.M|) ./^^ (60) 
Finally, for the box diagrams, we encounter rank-one and rank-two tensor integrals which 
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FIG. 8: Percentage difference between tlie Improved Born Approximation and the exact one-loop 
SQCD calculation for pp — )• 6/i as a function of mn at y/s = 7 TeV and tan f3 = 40. 



are written in terms of the Passarino-Vehmann coefficients as: 



p'tDu + p'^Du + p'sDu] ■ (61) 



167r^ 



D'-'iplplplpl {Pi+P2)\ {p2+P,f;MlMlMl) ^ 



I 



IGvr 



g^^ Dqq + tensor structures not needed here > . (62) 



Appendix B: One-Loop Results 

In this appendix we give the non-zero contributions of the individual diagrams in terms 
of the basis functions of Eq. [2D] and the decompositions of Eq. [22J The contributions 
proportional to m;, tan /3 are new and were not included in the results of Ref . [23i] . Although 
we specialize to the case of the lightest Higgs boson, h, our results are easily generalized 
to the heavier neutral Higgs boson, H, and so the Feynman diagrams in this appendix are 
shown for (pi = h,H. 
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FIG. 9: Total cross section for pp — )• b{b)h production including NLO QCD and SQCD corrections 
(dotted lines) as a function of renormalization/factorization scale using CTEQ6m (black) and 
MSTW2008 NLO (red) PDFs. We take Mg = 1 TeV and the remaining MSSM parameters as in 

Fig. m 

The self-energy diagrams of Fig. [TTJ 

i=i ^ ^ 
- -lt(-lY^^B„{Ml) (63) 

1=1 

where we have have used the shorthand notation for the arguments of Passarino-Veltman 
functions, B„ i (m?) = B„ i («;M|,M?). 

i=l ^ 

xg) = -lt(-ir^So(M^) (64) 

1=1 

andSo,! (m^.)-5o,i (^;M|,M^) 
The vertex functions of Fig. [T2j 
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FIG. 10: Total cross section for pp — t- h{h)h production including NLO QCD and SQCD corrections 
as a function of the factorization scale using MSTW2008 NLO PDFs. We take Mg = I TeV and 
the remaining MSSM parameters as in Fig. SI 




b cf b \ 



FIG. 11: Self-energy diagrams, Si and 52. 



Diagram Vi: 



X 



» 

Vi 



i=l 



t 



. 2 

^v! = -qJ2{^ (^12 + C23) + 2^24 - (-1)* 2mhMgS,-, {Co + Cu)} [m} 



i=l 
2 



X 



(1) 

Vi 



2 ^ 1(^12 + (^23 



(65) 



22 




FIG. 12: Virtual diagrams, Vi and V2- 



where Co,ii,i2,23,24 (m^J = Co,ii,i2,23,24 (O, 0, t; M?, M^^, M^-^ . 
Diagram V2: 

i=l 

4? = ^ E (-1)^ ^^^5^26 (^0 + Cu) (m^J (66) 

i 

where Co,ii,i2,23,24 (m^J = ^0,11,12,23,24 (o, 0, M|, M? , M? ) . 
The vertex functions of Fig. [T31 




FIG. 13: Virtual diagrams, V3 and V4. 
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Diagram V3: 

4? = y t {Cu + C.3 - (-!)■ (Co + C.)| (m|) 

i=l ^ 

= -2 5^ {^^^° - 2 (1 - e) C24 - (-1)^ 2m,M^S2^Ci2} (m?) 

4f = -3«E + C,, - (-1)^ ^^^^^ (Co + cA {mi) 
i=i ^ ^ 

i=l 
2 

4? = -3$^(-irm,M^S25{Co + Ci2}(M?) 

where Co,ii,i2,23,24 (m^J = ^0,11,12,23,24 (o, 0, t; M|, M|, M?) . 
Diagram V4: 

o 2 

4? = ~2 5^ {^^^° - 2 (1 - e) C24 - s (C12 + C23) + (-1)^ 2m,MgS^j^C^] 

i=l 

<^ = -3uE + - (-1)^ ^^^^ + C7i2)| (m|) 
j=i ^ ■' 

4? = -^E(-ir^.M,.,,Co(M|) 

i=l 

2 

4? = 3$^(-irm,M,.2^{Co + Ci2}(M?) 

i=l 

where Co,ii,i2,23,24 (m^J = Co,ii,i2,23,24 (0, 0, s- Ml M|, M? ) . 
The vertex functions of Fig. [Til 




FIG. 14: Virtual diagrams, V5 and Vg. 
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Diagram V5: 

2 

4? = ^^^> E C'/...5.,C^i2 (m|,M^J (69) 

where ^0,11,12,23,24 (m|^, = ^0,11,12,23,24 (o, M^, t; M|, M|^, M^) , the squark mixing 

matrix is defined, 



(70) 



«21 «22 y y C2b ^^2b 

and the hght Higgs-squark-squark couphngs Ch,ij, are normahzed with respect to the Higgs- 
quark-quark couphng[2|, 

Cmi + C^..22 = 4m, + ^-Mlifj^ (71) 
C^.,11 - C^.,22 = 2Y,s,i + ^c,^ (4 - 2Q,sl.) (72) 
C,,u = C^.,21 = Y,c,, - ^s,, {4 - 2Q^4) , (73) 



= sin = 1 - M^/Ml and is defined below Eq. E 
Diagram Vg: 

4 



= 3 5Z C7Mj{<^..^bC'ii + «i.^§C'o}(Mf,Mf^. 



< = =XIS = (74) 

where C„,„,,2,23,24 = Co,n,i2,23,24 (o,M;,s;Af|,M|_,M?) . 

The box diagram of Fig. [T3 



< = ^ E a.-^/.,.{^o + A3}(M|,Mf^, 

M=l,2 

<^ = E a,C,,,D,,{Ml,Ml) 

X^^l = 3M~gU E a,,CH.j{Dn-Di^s}[Ml,Ml 

i,i=l,2 

<^ = E ^^^C,,., {M|Do - 2Doo} (m^, M|) (75) 



«,i=i,2 
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FIG. 15: Box diagram, Bi. 



where, Do [Ml Ml ) ^ Do ( 0, 0, 0, Ml s, t; Ml Ml Ml M? 



The box diagram of Fig. [121 



nrsiTBTSTsy 




FIG. 16: Box diagram, i?2- 



Diagram B2: 

= E a,C,,,ADo + Dn}{MlMl) 

< = ^ E a,C,,,{Do + Dn}{MlMl) 

XS = ^I: a.,C,,.,{Du-D,,}[MlMl) 

= E ^^.-^^-^oo (m^,MJ) (76) 

where Do (^f^, ) = Do (o, 0, 0, Ml u, s; Ml Ml M?,Ml^y 
The box diagram of Fig. [T71 
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b 



b 




b ff« <l>i 

FIG. 17: Box diagram, B^. 



Diagram B^: 



X 



it) 



X 



(2) 
^3 



3 



(77) 



where Do (m^, J = Do (o, 0, 0, Ml u, t; M?, M?, M?,M?^^. 

The vertex and external wavefunction counter terms, Eq. [29l along with the subtraction 
of Eq. [32l give the counterterm of Eq. 

Air 



Snib 
rUb 

2 



i=l 



2MgYbI{M-,^,M-,^,Mg) + Y(- (-1)' ^^bS,iB'o + 2mlB[ ) (0; M|, M^^, 



(78) 



Note that the counterterm contains no large tan f3 enhanced contribution. 



Appendix C: Definitions 



In this appendix we define the f unctions used in the expansions of the Passarino-Veltman 



integrals in the maximum and minimum mixing scenarios, where R = ^ in the maximal 



Ms 
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mixing scenario, and Ri 

MR) 
h{R) 

h{R) 
h{R) 

hi (i?i,i?2,n) 



Mr 



Further, 



in the minimal mixing scenario: 



{1-R^y 
3 

(1 - R-^f 
4 

(1 - R^Y 
5 



[l-R'^ + R^ogR^] 
[l-R^ + 2RHogR^] 



(1 - R^f 

Rl Y^ogRl 



1 + ^R^- 3R^ + \r^ + log R^ 
^ - 4i?=^ + 4i?^ - \r^ - log R^ 



\-R\) \-R\ 
.^■ + 2 



Rl Y logi?i 



1 — R2 J 1 — R2 



/i2 (-Rl, R2) — 



-E(-i) 

Rl + Rl-2 



+ 



(1 - (1 - Rl-R 



2^2 



log Ri 



R^ + — 2R^ 
(1-^1) 



^logi?^ 



d/.(a;) 














1 d/,( 





dx2 



(79) 



(80) 



a;=l/i? 
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